Holographic Symmetries and Generalized Order Parameters for Topological Matter 
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We introduce a universally applicable method, based on the bond-algebraic theory ol dualities, 
to search for generalized order parameters in disparate systems including non-Landau systems with 
topological order. A key notion that we advance is that ol holographic symmetry. It reflects situations 
wherein global symmetries become, under a duality mapping, symmetries that act solely on the 
system's boundary. Holographic symmetries are naturally related to edge modes and localization. 
The utility of our approach is illustrated by systematically deriving generalized order parameters 
for pure and matter-coupled Abelian gauge theories, and for some models of topological matter. 
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Introduction. — Landau's concept of an order parame- 
ter (OP) and spontaneous symmetry breaking are central 
in physics [1]. In systems with long-range Landau orders, 
two-point correlation functions of an OP field 0(r), in 
their large distance limit, tend to a finite (i.e., non-zero) 

value, lim| I ._ r /|^ 00 lim A rd^oo(C( r )C^( r ')) 0> where N 
is the linear size of the d-dimensional system, and 0(r) 
is local in the (spatial) variable r. It is in Landau's spirit 
to use the OP as a macroscopic variable characterizing 
the ordered phase and as an indicator of a possible phase 
transition (classical or quantum) to a disordered state 
where the OP becomes zero. 

There is much experience, including systematic meth- 
ods [2, 3], for deriving Landau OPs and their effective 
field theories [1]. Landau's ideas of a (local) OP cannot 
be extended to topological states of matter as, by defini- 
tion [4, 5], these lie beyond Landau's paradigm. However, 
the notion of long-range order, or the design of a witness 
correlator (i.e., a correlator discerning the existence of 
various phases and related transitions), can be extended 
to topological phases - phases that can only be mean- 
ingfully examined by non-local probes [5]. Topological 
orders appear in gauge theories, quantum Hall and spin 
liquid states (when defined as deconfmed phases of emer- 
gent gauge theories [6]), including well-studied exactly 
solvable models [7, 8]. 

In this paper we demonstrate that generalized non- 
local OPs may diagnose topological phases of matter. 
Most importantly, we outline a method based on bond- 
algebraic duality mappings to search systematically for 
generalized OPs. Dualities have the striking capability of 
mapping Landau to topological orders and viceversa for 
essentially two reasons. First, dualities in general rep- 
resent non-local transformations of elementary degrees 
of freedom [9] and may even perform transmutation of 
statistics [10]. Second, bond-algebra techniques [10-12] 
allow for the generation of dualities in finite and infinite 
size systems. As we will show, in systems with a bound- 
ary, dualities realize a form of holography [13] capable 



of transforming a global symmetry that may drive spon- 
taneous symmetry breakdown into a boundary symme- 
try. We term these distinguished boundary symmetries 
holographic. They are under suitable further conditions 
connected to edge (boundary) states. To illustrate the 
method, we derive explicitly a (non-local) witness corre- 
lator and a generalized OP, suited to diagnose the transi- 
tion between deconfmed and confined phases of matter- 
coupled gauge theories, undetectable by standard OPs or 
Wilson loops. Other examples are reported in Ref [18]. 

The search for generalized order parameters. — A nat- 
ural mathematical language to describe a physical system 
is that for which the system's degrees of freedom couple 
locally. This simple observation is key to understanding 
that topological order is a property of a state(s) relative 
to the algebra of observables (defining the language) used 
to probe the system experimentally [5]. In the language 
in which the system is topologically ordered, it is also 
robust (at zero temperature [16]) against perturbations 
local in that language. Spectral properties are invari- 
ant under unitary transformations of the local Hamilto- 
nian H governing the system, H i-> UHU^. If UHU' 
corresponds to a sensible local theory then the unitary 
transformation U establishes a duality [10]. A duality 
may map a system that displays topological order to one 
that does not [5]. Dualities for several of Kitacv's models 
[7, 8, 17] epitomize this idea [5, 12, 16]. 

Since dualities are unitary transformations (or, more 
generally, partial isometries) [10] they cannot in general 
change a phase diagram, only its interpretation. This 
leads to a central point of our work: A duality mapping 
a Landau to a topologically ordered system must map the 
Landau OP to a generalized OP characterizing the topo- 
logical order. Our method for searching for generalized 
OPs, combines this observation with the advantages of 
the bond-algebraic theory of dualities [10]. In this frame- 
work, dualities in arbitrary size (finite or infinite) systems 
can be systematically searched for as alternative local 
representations of bond algebras of interactions associated 
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to a Hamiltonian H . Hence it is possible for any system 
possessing topological order to systematically search for 
a duality mapping it to a Landau order. When a dual 
Landau theory is found, the dual system's OP can be 
mapped back to obtain a generalized OP for the topo- 
logically ordered system. In what follows and in Ref. 
[18], we study various quantum gauge and topologically 
ordered theories, and their duals, to illustrate our ideas. 

Holographic symmetries and edge states: the gauged 
Kitaev wire. — We next illustrate the concept of holo- 
graphic symmetry and its relation to generalized OPs 
and edge modes. Consider the Kitaev wire Hamiltonian 
[17] with open boundary conditions, here generalized to 
include a Z2 gauge field (termed the gauged Kitaev wire) , 
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where a m = <4ri> b m = 6^ denote two Major ana fcrmions 
({a m ,a n } = 2S mn = {b m ,b n }, {a m ,b n } = 0) placed on 
each site of an open chain with N sites. The Pauli matri- 
ces a? m .iy a = x,z, placed on the links (m; 1) connecting 
sites m and m + 1 represent a Z2 gauge field. For the 
gauged Kitaev wire, fermionic parity is obtained as the 
product of the local (gauge) Z2 symmetries ibiaia^.^, 
a (N-i;i) ib NaN, and o-f m . 1) ^m+ia m+ iCT^ m+1;1) (m = 
1, ■ • • ,N — 2). Just like the standard Kitaev wire, 
HQK[h = 0] has two free edge modes a± and b^. 

The gauged Kitaev wire holds two important dualities. 
It is dual to the one-dimensional Z2 Higgs model [19] 
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with Pauli matrices erf placed on sites i. Moreover, the 
gauge-reducing [10] duality mapping $d 
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defined on (N + 1) sites. The fermionic parity P maps 
to a holographic symmetry under this duality, since 

P = T\m=i ib ™ a ™ a i a N+v i- e -' tne product of 

two (commuting) boundary symmetries. Holography is 
a relational phenomenon (see [18]). A duality that un- 
covers a holographic symmetry links a global (higher- 
dimensional) symmetry of a system to a boundary (lower- 
dimensional) symmetry of its dual. Boundary symme- 
tries need not in general be duals of global symmetries. 



What is the physical consequence of having an holo- 
graphic symmetry? Consider the not uncommon situa- 
tion in which the holographic symmetry is supplemented 
by an additional (non-commuting) boundary symmetry 
in some region of the phase diagram. By definition, holo- 
graphic symmetries are boundary symmetries which are 
dual to global symmetries. Thus, global symmetries link- 
ing degenerate states (and properties in the broken sym- 
metry phase) in the dual system have imprints in their 
holographic counterparts. Then, the many-body level de- 
generacy of the ground state may be ascribed to bound- 
ary effects. If the couplings are now changed, the ground 
state degeneracy may get removed, together with some 
boundary symmetries. However, so long as the system re- 
mains in a topological phase dual to the (broken symme- 
try) ordered phase, the low energy state splitting will be 
exponentially small in the system size, so that in the ther- 
modynamic limit ground state degeneracy is restored. 

The language providing the most local operator de- 
scription of the ground-state manifold is the one realiz- 
ing the edge modes, which are expected to be exponen- 
tially localized to the boundary. Thus, as long as the 
thermodynamic-limit degeneracy remains, a suitable lo- 
cal probe will detect localization on the boundary for 
those states. Conversely, non-commuting edge mode op- 
erators in a gapped phase reflect the existence of low- 
energy many-body states with energy splittings vanish- 
ing exponentially in the system size. Many-body (zero- 
energy) edge states are thus simply a natural consequence 
of a degenerate ground state manifold in a gapped sys- 
tem. They are witnesses of an ordered (degenerate) phase 
described in a most local language. Note that boundary 
operators that commute with the Hamiltonian at special 
values of the coupling(s) are a necessary but not sufficient 
condition to realize exact (zero-energy) edge modes. 

The duality — >■ -Z?gk maps a global symmetry of 
Hh[k, h = 0] to a boundary symmetry of Hq^[k, h = 0], 
i.e., of ■ - - crfj_ 1 cr^ ~~ ^ 1>n, and one boundary symmetry 
to another, erf —> a%. If we now turn on h < J, keeping 
k = 0, the edge mode operators ai,6j\r, evolve respec- 
tively into T 1 = ^(-W^dd 1 ^!)) and 

T 2 = ELi(V^)^ m t m (nf;>( s; i))- These modes 
are exponentially localized as long as the system is in 
the ordered gapped phase within a gauge sector [23]. 
The Majorana language affords a local boundary de- 
scription of these (partly non-local in the Higgs lan- 
guage) zero-energy modes. For h > J, and/or k > 0, 
the ground state is unique, even in the thermodynamic 
limit, as we learn from the phase diagram of the one- 
dimensional Higgs model [19]. Hence the zero-energy 
modes disappear together with the ground-state degen- 
eracy. For k > 0, they disappear despite the fact that 
fermionic parity remains an exact symmetry and can- 
not be spontaneously broken [21]. Consider now Hq K of 
Eq. (4). At ft, = 0, it has zero-energy edge mode oper- 
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ators af,af,a^ +1 ,a^ +1 . For h > 0, and k = 0, two of 
these remain unchanged, and the other two evolve into 
Si = of + ZZZl(h/J) m vl(n7=2°s)< + i and S 2 = 

4 +1 + EL J (V^^aiL+i o f)4 + i- These be - 

have just as their Majorana relatives, yet they are recog- 
nized as non-local. The Majorana language distinguishes 
itself as the most local one for zero modes. 

To obtain a generalized OP for the gauged Ki- 
taev wire, notice that Hq K [k = 0] reduces to the 
transverse-field Ising (Tl) chain. Hence it exhibits 
a second-order phase transition at J = h, K = 0. 
For Hq K [k = 0], this transition is witnessed by the 
Landau OP correlator lim^_ J |_ ! . 00 {Tl | crfcrj 1 |TI). (From 
now on | label) represents the ground state of iT| a bel)- 
Our duality maps this correlator back to a generalized 
OP for the gauged Kitaev wire, the string correlator 
lim| i _ J -|^ 00 (GK|i6 i a l z6 l+ iaj + i • • -ibjaj\GK). 

Generalized OPs in higher- dimensional theories — We 
next show how to systematically derive generalized OPs 
in higher space dimensions. Our main goal is to illus- 
trate the methodology in the challenging case of the 
Abclian (C/(l)) matter-coupled gauge (Higgs) theory. 
Previous works [14, 15] conjectured generalized OPs for 
matter-coupled gauge theories and were numerically im- 
plemented, for instance, in Ref. [6]. Unfortunately, a 
systematic mathematical derivation was missing and this 
is what our work is about. Our (non-local) witness cor- 
relator for the Higgs model turns out to be the one con- 
jectured in Ref. [14]. In Ref. [18], we study several other 
examples (displaying also holographic symmetries), in- 
cluding Ising and Z p gauge and Higgs theories, the Z p 
extended toric code [20] as an interesting example of 
topological order, and the XY model on the frustrated 
Kagome lattice. Non-Abelian extensions of our ideas 
based on Ref. [11] are currently under investigation. 
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FIG. 1. The Z gauge theory exactly dual to the quantum XY 
model must satisfy special boundary conditions and possesses 
a boundary symmetry. The lattice corresponding to the XY 
model is shown in thick lines, for N = 4. 

To derive the generalized OP for the Abclian Higgs the- 
ory, our starting point is the XY model defined in terms 



of continuous U(l) degrees of freedom ,s r = e~ ,er , 9 r e 
[0, 27r), placed at sites r = ie\ + je 2 = of a square 
lattice. The model's Hamiltonian reads 

N 

H XY = hJ2 (5) 

j N-l N N N-l 

+ 2 [ E E + EE S ^ra) + h.c] , 

i=i j=i i=i j=i 

with L r = —id/d8 r , and S^.^ = s r s^ +e ^. The XY 
model is dual to a Z (solid-on-solid like) gauge theory also 
defined on a square lattice, but with degrees of freedom 
X and R associated to links (r; /i = 1,2). (In matter- 
coupled gauge theories we will also have operators act- 
ing on sites r.) These operators satisfy X\m) = m\m), 
R\m) = \m — 1), R' \m) = \m + 1), with raeZ, and com- 
mute on different links (and/or sites). Then, the exact 
dual of Hxy for open boundary conditions reads 

N 

Hzg =hJ2 b 2 {id) (6) 
i,j=i 

j N N N N 

+ 2 [ E E R ^m + E E + ■ 

i=2 j=l i=l j=2 

We will call system indices the link indices (i, j; fi = 1, 2) 
labeling R operators that explicitly appear in Hzg, and 
extra indices the remaining link indices. In the bulk, the 
plaquette operator b^j-) reads 

b (i,j) = X (*,ja) + X (i+l,j;2) ~ X (i,j + l;l) ~ X (i,j;2)- (?) 

On the lattice boundary, the plaquette operators are set 
by two rules: (i) b {1 ^ N) = X( ljN . 1 )-X ( 2 iN .q-X( liN+1 . 1 y 
Thus, n) involves one degree of freedom Xn N+ i ; i) la- 
belled by an extra link index, (ii) The remaining bound- 
ary plaquettes are determined by Eq. (7) provided opera- 
tors labelled by extra link indices are omitted. With these 
definitions in tow, the mapping of bonds 

b {itj) A L(ij) , l<i,j<N, (8) 

R d,r,i) ^ S h-v,2 V l<i<N,2<j<N, 

R(i,j;2) ^ 2 < * < N > 1 < 3 < N, 

implements the duality transformation Hzg H\\. 
As the operators R^ jv+i-i) > -^n at+i-i) ^° no ^ a PP ear m 
Hzg, the operator Y( lv /v+i ; i) constitutes a boundary sym- 
metry of Hzg- Similar to the duality between the one- 
dimensional theories of Eqs. (1) and (4), this is a gauge- 
reducing duality. The gauge symmetries of Hzg > given by 

= R ^,j-l) R (i,j-,2) R \i-l,j-l) R \i, -l-2)^ 2 - j - N > 

are removed by the mapping since Au j\ 1. 

In the thermodynamic (N — > oo) limit, the strongly- 
coupled ( J 3> h) phase of the XY model displays sponta- 
neous symmetry breakdown of its global U(l) symmetry 
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with generator L\y = J2i j—i as evinced by a non- 

vanishing (XY|s r sJ,, |XY) in the limit |r — r'| — > oo. By 
virtue of being dual to the XY system, the gauge theory 
displays a non-analyticity in its ground state energy as h 
is varied and its symmetry is broken. However, the phase 
transition in the gauge theory cannot be characterized 
by a local OP. So, how can the duality connecting the 
two models bridge the drastic gap separating the physi- 
cal interpretation of their common phase diagram? The 
answer lies in our notion of holography, since 

N 
i,j=l 

Thus, the global symmetry of the XY model is holograph- 
ically dual to the (local) boundary symmetry Xn t N+i-2) 
of its dual gauge theory and cannot be spontaneously bro- 
ken in this dual theory [21]. This is how holographic sym- 
metries explain the non-Landau nature of critical tran- 
sitions in the Z gauge theory. There are no edge modes 
nor localization associated with this holographic symme- 
try as the ordered phase of the XY model is gapless. 

We now derive a generalized OP for the Z gauge theory. 
Let r be an oriented path from r to r' made of directed 
links ( € T, and we adopt the convention that Si = S7 r;(tt ) 
if I points from r to r + e M , or Si = S^ r \ if I points 

oppositely from r + e^ to r. Then s r .s^, = Yiier Also 
let T* denote the set of links I* such that &d(Ri*) = Si 
(r* need not be continuous, see Fig. 2). Then 

(ZG| J] R t .\ZG) A (XY| Sr 4,|XY>, (10) 

and so the string correlator on the left-hand side is a gen- 
eralized OP for the Z gauge theory, displaying long-range 
order in the ordered phase. On a closed path, Oj*er* Rl * 
reduces to a product of gauge symmetries. 
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FIG. 2. Dual sets of links F* and V. 

Finally, we couple the Z gauge theory to a Z matter 
field (defined on sites r), Hzh = Hzq + Hm, with 

H M =J2 [ X ( R r + R t) + K l lr;ri] » (11) 

r il=l,2 

and l(r;n) = X r+e ^ — qXi r . t)1 \ — X r . The resulting matter- 
coupled theory Hzh is dual to the Abelian Higgs model 



[19] with Hamiltonian 

H AH =^[x(B r + Bl)+hL 2 r (12) 

r 

+ EK,+^,+ s S)))'- 

Here S$j} •> = s r s ( r /i) s r+e M includes a coupling with in- 
teger charge q to the U(l) gauge field = e~ tB ( T ^ , 

S (r;/i) = e ' 9S(riM) i and B r = S( I .. 1 )S( I . +ei . 2 ) s ( r . + e 2 ;l) S (j-;2)- 

The correspondence between the two models, established 
by the mapping of bonds 
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which holds only on physical gauge-invariant states. The 
reason is that $d preserves all commutation relations 
while "trivializing" all gauge symmetries. More pre- 
cisely, Hzh's gauge symmetries G r = i? r ^4 r map to 
$d(G r ) = 1, while -//ah's gauge generators g,. = + 

L (r;2) - ^(r-ei;l) - £(r-e 2 ;2) - map to ^(flv) = 

as follows from Eqs. (13) ("l*^ 1 is the mapping obtained 
from Eqs. (13) by reversing all the arrows). 

If the Z matter field is weakly coupled to the Z gauge 
field, the string correlator of Eq. (10) will still change an- 
alytic behavior across transitions. Then from Eqs. (13) 

(ZH| J] i?HZH) A (AH| Sr 4,I]>/|AH), (14) 

we obtain a witness correlator for the Abelian Higgs 
model that reduces to a Wilson loop on closed contours 
(r = r') (here s^ = % if a link I points from r to 

r + e M and = s?j[ % otherwise). This non-local correla- 
tor is directly related to intuitively motivated generalized 
OPs like (AH| Sr 4, n JG r4|AH)/(AH|n iG r c *?|AH> con- 
jectured in earlier work [6, 14, 15] (Tc denotes a closed 
loop roughly twice as long as T and containing it). 

Outlook. — As demonstrated, holographic symmetries 
and generalized OPs appear in numerous systems once 
boundary conditions are properly accounted for in the 
framework of bond-algebraic dualities. By providing a 
systematic methodology and many examples, our results 
might bring the theory of generalized OPs and topolog- 
ical orders to a new level of development closer to that 
of Landau's theory. More key problems need to be tack- 
led. First, the sufficient conditions under which a given 
topological order may be mapped to a Landau order and 
viceversa should be understood. Second, the problem 
of associating effective field theories to generalized OPs 
should be studied systematically. 
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SUPPLEMENTAL MATERIAL 

Introduction. — In the main text, we illustrated key 
ideas by examining several examples: (i) the quantum 
one-dimensional Ising matter coupled gauge theory (the 
Z2 Higgs model) and its dual quantum spin chain in an 
external field, (ii) a two-dimensional quantum XY model 
and its solid-on-solid (Z) type gauge theory dual, and 
(iii) the quantum U(l) matter coupled (or Abelian Higgs) 
gauge theory and its Z matter field coupled gauge dual. 

In what follows, we illustrate how our concepts can 
be similarly worked out in other examples: (iv) the two- 
dimensional transverse field Ising model and its quantum 
Ising gauge theory dual, (v) the self-dual two-dimensional 
quantum Z2 Higgs model, (vi) the two-dimensional quan- 
tum p-clock model and its Z p gauge theory dual, (vii) the 
two-dimensional Z p Higgs theory which we earlier illus- 
trated to be self-dual for all p, (viii) the two-dimensional 
X p Higgs theory and its Extended Toric Code dual, and 
(ix) the quantum XY model on the kagome lattice and 
its Z gauge theory dual. Taken together, these examples 
demonstrate how holographic symmetries and general- 
ized OPs and associated correlators can be systematically 
derived in disparate theories. 

The two-dimensional transverse field Ising model and 
its quantum Ising gauge theory dual. — As in the main 
text, we will perform dualities for finite size systems and 
explicitly label the Cartesian coordinates 1 < i, j < N of 
sites r — (i,j) on an N x N square lattice. The planar 
transverse-field Ising model for such any N x N square 
lattice system is given by 

N 

ff, = -h (15) 

«')j=i 

JV-l N N JV-1 

i—1 j — 1 i—1 j — 1 

The first term denotes the transverse field along the x 
direction while last two terms correspond to Ising inter- 
actions (along the internal o~ z direction) associated with 
horizontal and vertical links of the lattice respectively. 
This transverse field Ising model is dual to (the gauge- 
invariant sector of) the Z2 Ising gauge theory in the in- 
finite size limit [1]. For open boundary conditions, we 
obtain via a bond-algebraic duality the exact dual of H\ 
which reads 

N 

H\g = ~hJ2 b ^j) (16) 

N N N N 

i—2 j=l i=l j=2 

Following our earlier conventions, the link indices 
(i,j;fi) = (r;(j,), p, — 1,2, that appear explicitly in H\q 



labeling spins a x will be termed allo wed indices, with all 
other link indices being forbidden. This convention will 
be invoked to define the plaquette product -B(ij) which 
will omit forbidden links. The plaquette operators will 
be of a similar nature to those of Fig. 1 of the main text. 
In the "bulk" , the plaquette operator Buj\ reads 



B 



cr (i,j;l) cr (i+l,j;2)0'(i,j+l;l) '(i,i;2)- 



(17) 



At the topmost left site (see Fig. 1), £?(i,jv) = 
a (i N-i^k N-2) a (i iv+i-iV H cncc B(x t N) involves one spin 
°fi iV+i-i) labelled by a forbidden link index (see Fig. 1). 
Elsewhere on the boundary of the square lattice, the re- 
maining boundary plaquettes are determined by Eq. (17) 
provided spins labelled by forbidden link indices are omit- 
ted. 

We now examine to the bond algebra- the set of 
algebraic relations between the interaction terms in the 
Hamiltonian of Eqs. (15, 16). This will allow us to 
establish a duality mapping between the transverse field 
Ising model and the Ising gauge theory. A natural par- 
tition of the transverse field Ising Hamiltonian H\ of Eq. 
(15) is into three type of bonds: (a) the transverse fields 
{°fi,3)}> ( b ) horizontal Ising interactions W^j^fi+ij)}, 
and (c) vertical Ising interactions { a (i j) a (i j+i)}- The 
bond algebra of this system is exhausted by the following 
relations: 

(1) All bonds of type (b) and (c) commute amongst 
themselves and each other. 

(2) Single site bonds of type (a) anticommute with 
the two-site bonds of type (b) [or (c)] whenever any two 
such bonds share a common site on the lattice; this is 
a trivial consequence of the anticommutation relation 
K,<t*} = 0. 

(3) Whenever bonds of type (a) share no common 
lattice sites with bonds of type (b) or (c) then any such 
bonds commute with one another. That this is so is 
readily seen as [af, a*,] =0 whenever r ^ r' . 

(4) The square of any bond is the identity operator 



(«r = «>. 



\2 



!))■ 



(5) The product of any horizontal/vertical bonds of 
types (b) and type (c) around any closed plaquette is 
equal to the identity operator. That is, as (a r ) 2 = 1 for 
any site r, for all 2 < i, j < N , 



1 = [o 



J °"(t-l,j-l) J 

x [°(»-lJ-l) 'fi,j-l)][ CT (iJ-l) CT fiJ)] 



(18) 



Upon removing gauge redundancies, an identical set of 
relations appears for the bonds of the Ising gauge theory 
of Eq. (16). To underscore this, we partition Eq. (16) 
into three types of bonds as follows: (a') the plaquette 
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operators {B^ ^}, (b') gauge fields erf -. 2 on vertical and 
(c') gauge fields erf ^ on longitudinal links. That similar- 
ity of algebraic relations underlies the duality mapping 
between the theories. That the two theories of Eqs. (f 5, 
16) are indeed dual to each other on an N X N square 
lattice (with of an size arbitrary N) is readily seen by the 
mapping of the bonds (a,b,c) — > (a',b',c'), or 



(19) 

j = 2,--,N, (20) 
» = 2,...,JV > (21) 



B {i,j) 


— > (i 




gj, , 




— > (i 




gd, , 




— > (i 



(i-lj/tij)' 



once the local gauge symmetries of the Ising gauge theo- 
ries are removed and set to unity, 



.4 



>(t-i.rA) a h-h2) — ^ !• ( 22 ) 

With the constraint of Eq. (22) implemented, redun- 
dant gauge degrees of freedom arc eliminated- applying 
a gauge transformation of the form of Eq. (22) does 
not lead to a new state. As the mapping of Eqs. (19, 
22) preserves all of the algebraic relations amongst the 
bonds in both theories of Eqs. (15, 16), this mapping im- 
plements an duality H\q — U H\ for any N x N square 
lattice. Similar to the main text, we briefly elaborate 
on the change of number of degrees of freedom once Eq. 
(22) is implemented and how this leads to an identical 
number of degrees of freedom in both the transverse field 
Ising and Ising gauge theory. The counting of the number 
of degrees of freedom is identical to that of the U(l) sys- 
tem examined in the main text. We repeat it anew here 
for the Ising systems at hand. The Hamiltonian H\q de- 
scribes 2( A — l)iV + 1 spins and possesses (N — l) 2 gauge 
symmetries which are removed by the constraint of Eq. 
(22). This leaves a total of 2(A-1)A+1-(A-1) 2 = N' 2 
independent spins. This is precisely the number of spins 
in H\. 

Now, we turn to a quintessential feature seen by this 
duality- that of the holographic symmetry. As, by virtue 
of our definitions above of forbidden indices, a z t N+1 . 1 - j 
makes an appearance only in 5n jy) ( see Fig. 1) an d 
^ljv+i-i) a PP ears nowhere in the Hamiltonian of Eq. 
(16), the operator a z x N+1 .^ constitutes a boundary sym- 
metry of H\c. What is this boundary symmetry's origin? 
The answer is afforded by a simple calculation, 



(l,iV+l;2) 



n 



B, 



(hi) 



(23) 



Thus, the boundary symmetry o~ z x N+1 . 2 ) * s holographic. 
As seen by Eq. (23), this symmetry is a "hidden" dual of 
the global 1i symmetry of the Ising model. Its recogni- 
tion became possible with the specific choice of operators 
and use of bond algebras which suited to treat exactly for 
finite size systems and their boundaries. 



We conclude our examination of this duality and its 
consequences by deriving the OP of the Ising gauge the- 
ory. As we will establish now, this OP will be given by 
string products- products of operators along connected 
links. To write this OP, we need to set a few prelim- 
inaries. Towards this end, as in the main text, we let 
r denote a path (set of concatenated links) in the Ising 
model's lattice, starting at site r\ and ending at r 2 so 
that ri(r,/i)er a r a r+e IJ , = <T r 1 a r 2 - We furthermore let 
(r*, e* ) denote that link in the lattice of the Ising gauge 

From these 



theory such that *d(of r . ie 
definitions, it follows that 



(IGI 



n 



|IG) 



l«|l>- 



(r*,e^*)Gr* 



(24) 

In Eq. (24), |IG) and |l) denote, respectively, the ground 
states of the Ising gauge (Eq. (16)) and Transverse field 
Ising model (Eq. (15)). Thus, putting all of the pieces to- 
gether, the (infinite separation limit of the) string corre- 
lator £?r* is the generalized OP of the Ising gauge theory. 
This OP is the dual of the Landau order as it appears in 
Eq. (??). 

The self-dual Z2 Higgs model in two dimensions. — We 
next examine the self-dual quantum Ising matter coupled 
gauge theory on the square lattice. Similar to its one- 
dimensional variant of Eq. (1), the Hamiltonian of this 
system is given by 



(hB r + Act*) 



E E W, s 

r fj,= l,2 



(r;M) 



Kcr, 



T (r; f j,)- a r+e ll 



(25) 
(26) 



Similar to the main text, In the Z2 Higgs Hamiltonian 
of Eq. (25), <7% ' z are Pauli operators that reside on the 
lattice site r (which play the role of quantum matter 
fields) while the quantum gauge fields t^ z ^ are Pauli 
operators that lie on a square lattice link that connects 
the site r with a neighboring site site r + [i along one of 
the two ({j, = 1, 2) square lattice directions. Following all 
that we have stated thus far, the reader can verify that 
the correspondence of bonds 



B r 



T (r;i) 



V;2) 



r._ e2 7"(t-- 



°V-ei T (r- 



(27) 

e 2 ;2) a ri ( 2 8) 
-e i; l)^. (29) 

along with analogous relations for the matter fields [6] 
proves that i?Higgs is self-dual in its gauge-invariant sec- 
tor. It is straightforward to check that T* maps to the 
path r just as before under the self-duality mapping of 
Eqs. (28,29). Hence, 



(h\ n <j h ) 

(r;n)er« 



<H|« 



n < 

(r;At)er 



fr;/*)l H >' 

(30) 
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where \H) denotes the ground state of the Higgs Hamil- 
tonian of Eq. (25). The gauge-invariant correlator on 
the right-hand side of Eq. (30) is related to known string 
OPs for Higgs theories [2, 3]. The reader is referred to 
these works for details concerning asymptotic behaviors 
of string operators in the confined and deconfined phases. 
For completeness, we briefly reiterate that in these OPs 
[2, 3], the expectation value of Eq. (30) is divided by 
string [2] (or Wilson loop [3] products which are derived 
from r. This division enables the extraction of the non- 
trivial asymptotic large Y behavior of the string corre- 
lator of Eq. (30); This factor is, equivalently, related to 
the normalization of the state generated, in an imaginary 
time representation, by the application of the string prod- 
uct on the ground state [3, 4]. Other relations to effective 
line tension resulting from combined bare line (borne by 
matter coupling) and surface (generated by pure gauge 
fields) tensions are discussed in [4]. 

The Z p Gauge Model. — We may place the Z 2 theories 
of the previous sections and the U(l) theories in the main 
part of the paper at the ends of an infinite sequence of 
models with Z p symmetry, p = 2, 3, • • ■ , 00. 

A "clock model" spin is associated to the group Z p 
is described by two unitary p x p matrices U, V satisfy- 
ing U p = t = V p and VU = ujUV, with u = c i2lxjv a 
pth root of unity [5]. When p = 2, these correspond to 
the Pauli operators a x ' z . Thus the Hamiltonian for the 
quantum p-clock model can be written as 



ff pc = -[^EE C7 w) c/ (Wj) ( 31 ) 



N-l N 



1=1 j=l 



j N N-l ^ N 



h.c. 



1=1 j=\ 



Similar to the Z2 theories, the p-clock is dual to (the 
gauge-invariant sector of) the Z p gauge theory 



JV N 



1=2 j=l 



j N N , N 

n EE V ^rA) + 9 E B (^) 



1=1 3=2 



h.c. 



M=l 



Following, verbatim, the same convention as earlier (as 
also marked in Fig. 1)) yet now applied to these oper- 
ators, the plaquette operator -B(jj) for plaquettes in the 
"bulk" is given by 



B 



U, 



t 



(33) 



On the boundary, = ^(i,JV;i)^2,JV;2)^(i,JV+i;i)- 

The remaining boundary plaquettes are determined by 
Eq. (33) provided operators labelled by forbidden link 
indices (those extend beyond the system boundaries) are 



omitted, see Fig. 1)). The correspondence of bonds 



B 



(ij) 



(34) 



V (itj . 2 ) A U {i _ ld) u\^ y i = 2, • • • , N, (36) 



establishes the (gauge- reducing) duality H p q ^> H p c- 
The operator lily N+1 .^ is an holographic symmetry of 
H p Q, since 



^(l,JV+l;2) 



i,j=l 



•I'd, 



N 

n v Lr 

i,3 = l 



(37) 



The operator on the right-hand side of Eq. (37) is the 
global Z p symmetry of the quantum p-clock model of Eq. 
(31). As we see from the left-hand side of Eq. (37), this 
operator is dual to the boundary symmetry of the gauge 



theory, ^ ( t ljN+1;2) . 
holographic. 

As we have in earlier examples, we now systematically 
devise an order parameter for the gauge theory. This 
Abelian gauge theory exhibits, as most others, both con- 
fined and deconfined phases. The program for methodi- 
cally deriving a generalized OP diagnostic will replicate 
our calculations in earlier examples. First, we write 
down a Landau correlator which directly captures the 
long range order in the Landau type system (in this case, 
the clock model of Eq. (31)). Relying on the duality 
between the quantum clock model and the quantum Z p 
gauge theory, we then determine the corresponding cor- 
relator which serves as a generalized OP for the gauge 
theory of Eq. (32). 

The Landau theory of the quantum p-clock exhibits 
long range order as seen by the finite value of the two- 
point correlation function, in the limit of infinite separa- 
tion, 



Thus, this symmetry, l/J 1 



(1,JV+1;2)' 



(pqt^ffj, . |pC), ^(i-i'Y + U-]') 2 -> co, 

(38) 

where | pC) denotes the ground state of the p-clock model 
of Eq. (31). We now map this two-point correlator onto 
a generalized OP for the gauge theory by invoking the 
duality that we found. Towards this end, our steps will 
replicate, yet once again, those that we invoked in earlier 
examples. Let Y be an oriented path from r to r' , and 
! 6 T its links oriented accordingly. Furthermore, we 
employ the convention that Ui = U r Ul +e ^ if I points 
from r to r + e M , and Ui = U}.U r + e if I points oppositely 
from r + e M to r. With this conventions, 



(39) 



If we define Y* as the set of links mapping to the path Y 
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under duality and $d(Vj») = Ui, then we will find that 
G = (pG\ [] ^.|pG) A ( P C\U r ul,\pC), (40) 

where |pG) is the ground state of the gauge theory H p q. 
In the confined phase of the gauge theory, the generalized 
OP Q assumes a non-zero value in the limit of infinite T* . 
By contrast, in its deconfincd phase, the gauge theory of 
Eq. (32) has a vanishing value of Q as V* is made asymp- 
totically large. This result constitutes a generalization of 
our earlier Ising (i.e., p = 2) results. Thus, similar to 
our earlier examples, putting all of the pieces together, 
the (infinite separation limit of the) string correlator Q is 
a generalized OP for the Z p gauge theory distinguishing 
confinement from deconfincmcnt, and derived from sim- 
ple considerations of spontaneous symmetry breakdown, 
duality, and holographic symmetries. 

The self- dual 1 V Higgs model. — We now extend the 
earlier Ising results to the Z p system, review the self- 
duality of the Z p gauge theory when it is coupled to Z p 
matter fields, and derive new string operators for this 
theory. This matter coupled gauge theory is the quantum 
Z p Higgs theory. That is, we will augment the (pure) 
gauge theory H p q of Eq. (32) by terms that minimally 
couple the gauge fields on the lattice links U r ^ to matter 
fields U r and 77 r +e M at sites that form the endpoints of 
those links as well as by additional matter only terms V r , 



V r + - 



E 

(11=1,2 



UrU (r,fi) U r+ efJ 



+ h.c. 



(41) 

From the minimal coupling term, where q is the power 
of the gauge field on the link that is coupled to the 
matter field, similar to the 77(1) example in the main 
text, it is evident that q plays the role of a charge. 
As U q = uCp-9)t, this charge can assume the values 
q = 1, ■ ■ • , Int\p/2] (where Int[-] denotes the integer par 
of its argument). On augmenting the Hamiltonian of 
Eq. (32) by the terms in Eq. (41), the resulting model 
Hamiltonian H p n- the (charge q) Z p Higgs model- is self- 
dual in its gauge invariant sector [6]. Specifically, this 
self-duality is given by the mapping 











V(r,l) 








V(r,2) 




U r - 


- ei U 



5ft 



21,1)^' 



for the gauge fields, and the transformation 



TT Tjl TV 



H v, 



(r-e 2 ,2)' 



^t^(r,2)^ r + e 2 ~~ ^ ^(t-ei,X)' 



(42) 
(43) 
(44) 

(45) 
(46) 
(47) 



The string operator in Eq. (40), becomes, on invoking 
Eqs. (43, 44), 

<pH| [] V t .\pH) ^ (pHl^n^'IpH). (48) 



In Eq. (48), |pH) denotes the ground state of the Higgs 
Hamiltonian (given by (H p q + Hmatter)- the sum of the 
Hamiltonians in Eqs. (32, 41)). Similar to our earlier 
conventions, in Eq. (48) we set Ui = Uf p n] if I points 

from r to r + e M , and Ui = U 



otherwise. 



As in our discussions for the other examples of Higgs 
theories concerning string OPs and their normalization, 
in the infinite separation (i.e., T) limit, the gauge- 
invariant correlator on the right-hand side is related to a 
normalized generalized OP proposed in Ref. [2, 3]. Taken 
on a closed path, the right-hand side of Eq. (48) reduces 
to a Wilson loop. 

The Z p Extended Toric Code model. — The Toric code 
model [7] and its extensions, e.g., [6, 8, 9] have gained 
much interest in recent years. Their original motivation 
was to serve as a caricature of toy models in which some 
of the basic concepts of topological quantum computing 
can be explained [7]. Since the model and its variants 
have been investigated heavily large on their own right. 
The original Toric code model had an Ising symmetry. In 
the current section, we will further employ its 1 p variant. 

In the current context, our interest in this system de- 
rives from the fact that the Z p Higgs model of the pre- 
vious section is dual to a generalization of the extended 
Z p Toric Code model (pETC). This model is given by 



pETC 



— B r 
2 

h.c. . 



(49) 
(50) 



on 



the 



square 



2 ) ' (r-ei,l) " (r-ea,2) 

Toric Code is recovered for p 



lattice, with A r = 
The standard extended 
2, in which case q — 1 is 



the only possible value. The Z p toric code discussed in 
Ref. [8] is recovered for q = 1. 

A gauge-reducing duality connecting H p n to i? P ETC is 
given by 



v r ' ^ r ' 

UrU (r,l) U r+ ei 



-% u q 



(51) 



for the matter fields. 



r7 t 77 9t 77 _l 77 9t 

while leaving the gauge fields invariant , B r B r 

and V( r>f A V( rtfl ). This duality removes degrees of 

freedom that lie at lattice sites. While the Hamiltonian 
fipETC has no global symmetries, the Z p Higgs model has 
one global symmetry. Similar to our earlier discussions, 
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the duality H p h — ^ H p ejc maps a global symmetry of 
H p h to a holographic symmetry of H p etc- 

The generalized OPs introduced in the previous sec- 
tion for the Z p Higgs model can be mapped to (dual) 
generalized OPs for the (self-dual) extended Toric Code 
model. Employing Eqs. (51) on the correlator of Eq. 
(48), 



( P H| ] ; y,.| P H> -% ( P ETC| ] ; VHpETC), (52) 



and 



(pHlC/^n^lpH) -% (pETClH^lpETC), (53) 



ier 



where |pETC) is the ground state of the extended toric 
code model of Eq. (49) and, as in Eq. (48), pH) denotes 
the ground state of the Higgs model. As H p u is self-dual 
in its gauge-invariant sector, i/pE~rc is indeed also self- 
dual. The two string operators of Eqs. (52,53) transform 
into each other under this self-duality transformation. 

The quantum XY model on the kagome lattice. — We 
now examine a frustrated spin system with a continuous 
symmetry: the quantum XY model on the open kagome 
lattice- [see Fig. 3, we will denote this lattice by Ajf], 



H: 



XYK 



» 21 

r£A K 



Li 



E k < 

{r,r') 



(54) 



The XY rotor at any kagome lattice site r has an orienta- 
tion 9 r and L r = ig§-- This XY Hamiltonian is invariant 
under global rotations generated by the total angular mo- 
mentum L = ^2reA L r - A duality maps this model to 
a solid-on-solid like (Z) gauge theory (similar to the one 
discussed in the main text) which now appears on the 
dice lattice A^, 

#x d yk= E Wi+rI)+ E ^(Ec- 1 )^) 2 - 

ledp 



leAi 



peAi 



(55) 

Here, I are the oriented links of the dice lattice, p its 
plaquettes, and (—1)' = ±1, with (—1)' = —1 only if I 
viewed as a link in has opposite orientation to that 
assigned by I £ dp (dp = oriented boundary of p). Sim- 
ilar to the examples that we discussed earlier, the exact 
dual possesses incomplete plaquette at the boundaries; 
this is seen in the lower panel of the figure. Under the 
duality transformation, the generator of the global ro- 
tational symmetry of the XY model maps onto a holo- 
graphic symmetry, 



L 



E 

r€A K 



L p 



E 

peA E 



ledp 



(-l) l X t = X v 



(56) 



Here, Xy is a degree of freedom placed at the boundary 
link marked with an open circle. Any other Xi enters 





ATA" 



FIG. 3. The quantum, U(l) symmetric XY model on the 
kagome lattice is holographically dual to a solid-on-solid-like 
gauge theory on the dice lattice. 



the twice sum, with opposite sign. This is then an exam- 
ple of holography for a continuous symmetry in a highly 
frustrated system. 
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